Defect production in non-linear quench across a quantum critical point 
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We show that the defect density n, for a slow non-linear power-law quench with a rate r _1 and an 
exponent a > 0, which takes the system through a critical point characterized by correlation length 
and dynamical critical exponents v and z, scales as n ~ r -««"i/(«*«+i) [ n „ ( aa (<»-i)/°</ r )<'d/(*»'+i)] 
if the quench takes the system across the critical point at time t = [t = to 0], where p is a non- 
universal constant and d is the system dimension. These scaling laws constitute the first theoretical 
results for defect production in non-linear quenches across quantum critical points and reproduce 
their well-known counterpart for a linear quench (a = 1) as a special case. We supplement our 
results with numerical studies of well-known models and suggest experiments to test our theory. 

PACS numbers: 73.43.Nq, 05.70.Jk, 64.60.Ht, 75.10.Jm 
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Quantum phase transitions have been studied exten- 
sively for several years Such transitions are accom- 
panied by diverging length and time scales [l[ leading to 
the absence of adiabaticity close to the quantum criti- 
cal point (QCP). Thus a time evolution of a parameter 
X(t) = Ao|i/r| Q sign(i) in the system Hamiltonian, char- 
acterized by a rate 1/r and an exponent a > 0, which 
takes such a system across a QCP located at A = A c , 
leads to a failure of the system to follow the instan- 
taneous ground state and hence to the production of 
defects @, 0). All the previous studies on such sys- 
tems have been restricted to linear quenches (a = 1) 
@, i, H, 0, S, @, EU El El ■ For such quenches, the density 
of these defects n ~ T -vd/(zv+i) ^ w ]j ere v anc [ z are 

correlation length and the dynamical critical exponents 
characterizing the critical point, and d is the system di- 



13, LLJ]. On the experimental side, trapped ul- 



mension 

tracold atoms in optical lattices provide possibilities of re- 
alization of several model quantum spin systems and are 
particularlysuitable for studying their non-equilibrium 
dynamics 15, [lB|. Experimental studies of defect pro- 
duction have indeed been undertaken for a spin-one Bose 
condensate 17\. Although these experiments can easily 
investigate non-linear quenches, there have been no the- 
oretical studies of such quench phenomena so far. 

In this letter, we show that a slow non-linear power- 
law quench, as discussed above, through a QCP leads 
to a density of defects which scales either as n ~ 

T -aud/(azu+l) Qr ag n _ ( ag (a-iya/ T }ud/(zv+l) ( wner e 

g is a non-universal model dependent constant), depend- 
ing on whether the quench parameter A vanishes or stays 
finite at the critical point. Such a scaling law for the de- 
fect density generalizes its earlier known counterpart for 
linear quenches 13] and thus constitutes a significant ex- 
tension of our understanding of quench dynamics across 
a QCP. Our results, to the best of our knowledge, also 
constitute the first theoretical investigation of defect pro- 
duction due to non- linear power-law quenches. We sup- 
plement our theoretical results with numerical studies of 



the one-dimensional Ising and Kitaev models, and also 
suggest realistic experiments to test our theory. 

We begin our analysis with a study of a model Hamil- 
tonian in d dimensions of the form 

H = J2 V^(fc) H(k;t) Tjj(k), 
fc 

H(k;t) = [A(t) + fe(fc)]r 3 + A(fc)T + + A*(fc)r_, (1) 

where 74=1,2,3 are the usual Pauli matrices with t± — 
(ri ± ivz)/2, b(k) and A(fc) are model dependent func- 
tions, and tp(k) = (c'^fc), c' 2 (fc)) represents fermionic op- 
erators. Many of the model systems with a QCP charac- 
terized by v — z = 1, such as the Ising and XY models 
0,0] in d — 1 and the extended Kitaev models [l8l.fl9l.l20j 
in d = 2, can be mapped onto such a fermionic Hamil- 
tonian via standard Jordan- Wigner transformation. We 
first consider the case where the system passes through 
a gapless point at t = and k — ko. Note that in this 
case both b(k) and A(fc) must vanish at k = kg. In what 
follows, we shall also assume that |A(fc)| ~ \k — ko\ and 
b(k) ~ \k— ko\ Zl at the critical point, where z\ > 1 so that 
E ~ \k — ko\ and z = 1. In the rest of the analysis, we 
set h= 1, and scale t — > tXo, t — > tAo, A(fc) — > A(fc)/Ao, 
and b(k) -> b(k)/X . 

The dynamics of the such a system is governed by 
the Schrodinger equation given by id t ^p{k) = H{k; t)ip(k) 
which leads to the equation governing the time evolution 
of ci(jfe) = c ' 1 (fc)e i / tdt '[l t '/^r^gn(t')+6(fe)] j 



dt 2 



2i 



sign(t) + b(k) 



d_ 



A(fc)| 2 Cl (fc) = 0.(2) 



Now we scale t — > tT a ^ a+1 "> so that Eq. becomes 
d 2 



dt 2 



2i 



d_ 

Tit 



|t| a sign(t) + b(k)r— 
+ |A(fe)| 2 r^r) ci(A) = 0. (3) 



2 



Noting that the system was in the ground state 
(ci(fc), 02(h)) — (1,0) at the beginning of the quench at 
t = —00, we find, using Eq. ([3]), that the defect proba- 
bility p(k) must be given by 

p(k) = fim \a(k,t)\ 2 = / \b(k)T^TT; \A(k)\ 2 T^fr 

' (4) 

where / is a function whose analytical form is not known 
for a / 1. Nevertheless, we note that for a slow quench 
(large r), p(k) becomes appreciable only when the instan- 
taneous energy gap SE = 2[(A(i) + b(k)) 2 + |A(fc)| 2 ] 1/2 
becomes small at some point of time during the quench. 
Consequently, / must vanish when either of its argu- 
ments are large. Thus for a slow quench (large r), the 
defect density n = J BZ d d k/ d f \b(k)r^ ■ \ A(k)\ 2 r^ 
(where BZ denotes the Brillouin zone) receives its main 
contribution from values of / near k = k$ where both 
b(k) and A(fc) vanish and can be written as (extend- 
ing the range of momentum integration to 00) n ~ 



Jd d k/(2TT) d f 



\k-ko\ 



*+ 1 ;|/c-fc rT-+ 1 



Now scal- 



ing 



k^ (k- fc )i- Q/(Q+1) , we find that 



n = t a+1 



d d k 



(27T) 

d d k 



^f{\k\ Zl r a(1 - Zl)l(a+1) -Ak\) 



(5) 



where in arriving at the last line, we have used z\ > 1 
and t — > 00 . (If z\ = 1, the integral in the first line 
is independent of r, so the scaling argument still holds). 
Note that for a = 1, Eq. ([5]) reduces to its counterpart 
for a linear quench [1 31 ] . 

Next we generalize our results for a critical point with 
arbitrary v and z. To this end, we consider a generic time 
dependent Hamiltonian Hi[t] = Hi[X(t)], whose states 
are labeled by \k) and |0) denotes the ground state. If 
there is a second order phase transition, the basis states 
change continuously with time during this evolution and 
can be written as \ip{t)} = J2k a k(^)\k[^(t)}} an d the de- 
fect density can be obtained in terms of the coefficients 
aj:(t) as n = J2k^o \ a k^ ~^ °°)l 2 so that one gets [l]| 



d d k 



(27T) 



dX 



(6) 



where 6u}^(X) = k>z{X) — u>q(X) are the instantaneous ex- 
citation energies, and we have replaced the sum over k 
by a rf-dimcnsional momentum integral. We note, fol- 
lowing Ref. [T3L that near a critical point, Slu^(X) = 
AF(A/\k\ z ), where A is the energy gap, z is the dy- 
namical critical exponent and F(x) ~ 1/x for large 
x. Also, since the quench term vanishes at the critical 
point, A ~ |A| QZI/ for a non-linear quench, one can write 
Suj^(X) = \X\ az "F'(\X\ azi, /\k\ z ) where F'(x) - 1/x for 



large x. Further, one has (fe|^|0) = \k\- z G(A/\k\ z ) 
near a critical point where G(0) is a constant. This al- 
lows us to write (fc|^|0) = ^^-G' [X azu /\k\ z ) where 
G'(0) is a constant [HEDl]. Putting these in Eq. © and 
changing the integration variables to r\ = T au /( azu + 1 )\k\ 
and £ = l^-^^A, we find that 



£j — aud/(<xzv-\-X) 



(7) 



where C is a non-universal number independent of r. 

Next we focus on the case where the quench term does 
not vanish at the QCP for k = kg. We again consider 
the Hamiltonian H(k) in Eq. (TT]), but now assume that 
the critical point is reached at t = to 7^ 0. This renders 
our previous scaling argument invalid since A(fco) = 
but b(ko) 7^ 0. In this situation, \to/r\ = g 1 /", where 
g = \b(ko)\ is a non-universal model dependent constant, 
so that the energy gap 5E may vanish at the critical 
point for k = ko- We now note that the most important 
contribution to the defect production comes from times 
near £o and from momenta near kg. Therefore, we expand 
the diagonal terms in H(k) about t = to and k = ko to 
obtain 



ag 



(a-l)/a 



t-tr 



+ b'(Sk)\T 3 



+ A(k)T + + A*(k)r- t(,(k), (8) 



where b'(Sk) represents all terms in the expansion of b(k) 
about k = ko and we have neglected all terms R n — (a ~ 
n + l){a-n + 2)...(a-l)g^- n+1 ^ a \(t-to)/T\ n sign{t)/n\ 
for n > 1 in the expansion of t/r about to. We shall 
justify neglecting these higher order terms shortly. 

Eq. (J5]) describes a linear quench of the system with 
Teff(a) = T/(ag( a ~ 1 ^ a ). Hence one can use the well- 
known results of Landau-Zener dynamics (2l| to write 
an expression for the defect density n = J BZ d d k p(k) = 
J BZ d d kexp[—Tr\A(k)\ 2 T e g(a)]. For a slow quench, the 
contribution to n comes from k near ko so that 



Toff (a) 



-d/2 



(a-l)/a 



/r 



d/2 



(9) 



Note that for a = 1, we get back the familiar result 



-d/2 



as a special case and the dependence of n on 



n ~ r 

the non-universal constant g vanishes. Also, since the 
quench is effectively linear, one can use the results of 
Ref. 13| to find the scaling of the defect density when 
the critical point at t = to is characterized by arbitrary 
v and z. 



(ag^-^/r) 



vd/(zv+l) 



(10) 



Next we justify neglecting higher order terms R n . We 
note that the significant contribution to n comes at times 




a 



FIG. 1: Variation of the defect density n with the quench 
exponent a for representative values of r = 10 (black solid 
line), r = 15 (red dashed line) and r = 20 (blue dotted line). 
A polynomial fit of the form n = aa b yields exponents which 
are very close to the theoretical result 1/2 for all values of r. 
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FIG. 2: Plot of In(n) vs ln(-r) for the ID Kitaev model for 
a — 2 (black solid line), a — 4 (red dotted line), a — 6 (blue 
dashed line) and a — 8 (green dash-dotted line). The slopes 
of these lines agree reasonably with the predicted theoretical 
values —a/(a + 1) as shown in the table. 



t when the energy levels of H'(k; t) (Eq. ([5])) for a given 
k are close to each other: (t — t )/r ~ A(fc). Also, 
for a slow quench, the contribution to the defect den- 
sity is significant only when p(k) is significant, ie., when 
|A(/c)| 2 ~ l/r c ff(a). Using these arguments, it is easy 
to see that R n /R n -i = (a — n + l)g~ 1/,a (< — t )/(n,T) ~ 
(a — n + l)/(ni/r). Thus we find that all higher order 
terms R n >i, which were neglected in arriving at Eq. (|9j), 
are unimportant in the limit of slow quench (large r). 

The scaling relations for the defect density n given by 
Eqs. (0 and (fTT)|) represent the central results of this 
letter. For such power-law quenches, unlike their linear 
counterpart, n depends crucially on whether the quench 
term vanishes at the critical point. For quenches which 
do not vanish at the critical point, n scales with the same 
exponent as that of a linear quench, but is character- 
ized by a modified non-universal effective rate T e g(a). If, 
however, the quench term itself vanishes at the critical 
point, we find that n scales with a novel a dependent 
exponent avdj{azv + 1). For a = 1, T c g(a) = t and 
avd/{azv + 1) = vd/{zv + 1); hence both Eqs. ([7]) and 
(|10p reproduce the well-known defect production law for 
linear quenches as a special case [13|. We note that the 
scaling of n will show a cross-over between the expres- 
sions given in Eqs. ([7|) and (I10| near some value of t = tq 
which can be found by equating these two expressions; 
thisyieldsTo ~ |6(fco)| _zly_1 ^ Q - Fora > 1, the scaling law 
will thus be given by Eq. [7] (Eq. [TUJl for r < (>)r . We 
also note here that our results do not apply to quenches 
which take a system through a critical line 13, 22|. 

We now supplement these analytical results with nu- 
merical studies of well-known models. The first model 



that we choose for this purpose is the one-dimensional 
Ising model in a transverse field with the Hamiltonian 
Rising = -J(J2(ij) SfSj-90 Y,i s f) where J is the near- 
est neighbor coupling and go is the dimensionless trans- 
verse field. A standard Jordan- Wigner transformation 
[l[ then maps -ffising to a free fermionic Hamiltonian 

HLnJJ = Ek4l(9o - cos(fc))r 3 + sin(fc)rity*. Thus 
a time variation go(t) = | i/r | Q sign (£) takes the system 
through two critical points at to = t(— t) where the en- 
ergy gap vanishes at ko — 0(vr) so that go — 1(— 1) at 
these points. Thus the defect production around both 
these critical points have the same r e gf(a) = r/ot and we 
expect (Eq. ©) the defect density to go as n ~ -fa. 
for a fixed r. To confirm this expectation, we solve 
the time-dependent Schrodinger equation idtip(k, t) — 
H[ sin ij){k,€) and compute the defect probability pk and 
hence n for fixed r and for several representative values 
of a > 1 . These values of a and r are chosen so that we 
are in the regime where all R n >i can be safely neglected. 
The plot of n as a function of a for r = 10, 15, and 20 
is shown in Fig. [TJ A fit to these curves yields expo- 
nents of 0.506 ± 0.006 (r = 10), 0.504 ± 0.004 (r = 15), 
and 0.505 ± 0.002 (r = 20) which are indeed remark- 
ably close to the theoretical value 1/2 predicted by Eq. 
([9]). The systematic positive deviations in the exponents 
comes from the neglected terms R n >i- We note that the 
range of a for which such deviation remains small grows 
with r, as expected from our theoretical prediction. 

Next, we consider the one-dimensional Kitaev model 
[20L [23, which has the Hamiltonian Hk = 

Eleven (^fSf+i + J 2 Sf S?^) , where the sum ex- 
tends over even sites i on the disconnected chains of 
the underlying hexagonal lattice, and Si denotes the 



4 



spin at site i. Such a model can be realized as 
the J3 = limit of the well-known Kitaev model 
and can be mapped, via a standard Jordan- Wigner 
transformation [2Cj, [24J, onto the fcrmionic Hamilto- 
nian H' K = 2 ^[ {-J- sin(fc)r 3 + J+ cos(fc)r 2 ) ipk = 
2 Ek ^k H '( k )^k where i/>(k) = (ci(fc), c 2 (fc)) are 
fcrmionic operators, < k < tt extends over half the 
Brillouin zone, J± = J\ ± J 2 , and we have chosen 
the lattice spacing to be unity. Here the time vari- 
ation J-(t) = J|t/r| Q sign(i), keeping J + fixed, takes 
the system through a single critical point at t = and 
fc = kg = 7r/2 which has f = z = 1. The defect density, 
according to Eq. ([7]), is therefore expected to scale as 



-a/(a+l) 



To check this prediction, we numerically 



solve the Schrodingcr equation idtip(k) — H' K (k; t)ip(k, t) 
and compute the defect density n — J Q dk/ / np{k) as a 
function of the quench rate r for a with fixed J+/ J = 1. 
A plot of ln(n) as a function of ln(r) for different values 
of a is shown in Fig. [2j The slope of these lines, as can 
be seen from Fig. [2j changes from —0.67 towards —1 as 
a increases from 2 towards larger values. This behavior 
is consistent with the prediction of Eq. ([7]) . The slopes 
of these lines also show excellent agreement with Eq. ([7]) 
as shown in the inset of Fig. [2] 

Experimental verification of our results may be 
achieved in several possible ways. First, there has been 
a concrete proposal for the realization of the Kitaev 
model using an optical lattice In such a realiza- 

tion, all the couplings can be independently tuned us- 
ing separate microwave radiations. In the proposed ex- 
periment, one needs to keep J3 = and vary J 1( - 2 ) = 
J(l ± |t/r| Q sign(t))/2 so that J+ remains constant while 
J_ varies in time. The variation of the defect density, 
which in the experimental setup would correspond to the 
bosons being in the wrong spin state, would then show 
the theoretically predicted power-law behavior (Eq. (|7f). 
Secondly, a similar quench experiment can be carried out 
with spin one bosons in a magnetic field described by an 
effective Hamiltonian H cS = c 2 n (S) 2 + ciB 2 (S 2 } [13] 
where c 2 < and n is the boson density. Such a system 
undergoes a quantum phase transition from the ferro- 
magnetic to polar condensate at B* = \/|c 2 |no/ci. A 
quench of the magnetic field B 2 — B 2 \t/r\ a thus would 
lead to scaling of defect density with an effective rate 
Teff(a) = T J /(a</ a_1 ' ) / Q ), where g = |c 2 |no/ci. A mea- 
surement of the dependence of the defect density n on a 
should therefore serve as test of prediction of Eq. (fTT))) . 
Finally, spin gap dimer compounds such as BaCuSi 2 06 
are known to undergo a singlet-triplet quantum phase 
transition at B c ~ 23. 5T which is known to be very 
well described by the mean-field exponents z = 2 and 
v = 2/3 [25| . Thus a non-linear quench of the magnetic 
field through its critical value B = B c + _Bo|i/r|"sign(i) 
should lead to scaling of defects n ~ T -6a/(4a+3) m ^ = 3 
In the experiment, the defect density would correspond 
to residual singlets in the final state which can be com- 



puted by measuring the total magnetization of the system 
immediately after the quench. 

To conclude, we have obtained general scaling laws 
of the defect density for an arbitrary power-law quench 
through a QCP which reproduce their linear counterpart 
as a special case. We have verified our theoretical predic- 
tion by numerical simulation of model systems and have 
suggested several possible experiments to test our results. 
Our results have been recently used to find the optimal 
passage through a QCP [26j. 
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